Correlation and response in a driven dissipative model 
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We consider a simple dissipative system with spatial structure in contact with a heat bath. The 
system always exhibits correlations except in the cases of zero and maximal dissipation. We explicitly 
calculate the correlation function and the nonlocal response function of the system and show that 
they have the same spatial dependence. Finally, we examine heat transfer in the model, which 
agrees qualitatively with simulations of vibrated granular gases. 
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Driven dissipative systems (DDS) occur in many dif- 
ferent contexts, from collections of macroscopic particles 
to biological systems. Such systems are intrinsically out 
of equilibrium, and need the input of energy in order 
to remain functional. One prototype DDS is a granu- 
lar gas|3, a collection of inelastic grains which dissipate 
energy through collisions; these have been examined ex- 
tensively, both experimentally0,0|5 numerically 0, H, 0, 
and analytically[2|. Because of the difficulty in treating 
granular gasses analytically, stochastic mean-field mod- 
els have been studied. One such is the Maxwell modelQ, 
which assumes a velocity independent collision rate and 
no spatial structure; this facilitates analytical calculation 
of quantities such as the velocity distribution function 
and its moments. 

The introduction of spatial dependence complicates 
and enriches the behavior, and may lead to correlations: 
for example, actual granular gasses exhibit spatial clus- 
tering and velocity correlations because of the dissipa- 
tive collisions, as is seen in simulational studies |inl |. 
Williams and MacKintosh^^ have shown numerically 
that correlations exist in a ID driven dissipative gas pro- 
vided the restitution coefficient is different from one, and 
Soto, et oi[l^ have used the BBGKY hierarchy to study 
the appearance of velocity correlations in inelastic hard- 
sphere systems. Baldassarri, et al\3], and Ben-Naim and 
Krapivsky^^ consider a lattice variant of the Maxwell 
model which they solve in the freely cooling case (no driv- 
ing); the latter authors calculate the spatially dependent 
velocity correlations which exhibit gaussian decay with 
distance. 

In this Rapid Communication wc study a model of 
a driven system with spatial structure: the constituent 
"particles" are constrained to lie on a ID lattice with 
nearest-neighbor coupling. Our main goal will be to un- 
derstand the connection between the system's dissipa- 
tive nature and spatial correlations. The system is cou- 
pled to a heat reservoir at temperature T, and the model 
is chosen so that it has a well-defined equilibrium limit 
for certain values of the system parameters. The mean- 
field version of the model, which has no spatial structure, 
can be solved exactly (l^, in the sense that all the mo- 
ments of the energy distribution may be calculated. For 



the model of this Rapid Communication, we calculate 
the two-point correlations of the system analytically, and 
demonstrate that non-zero correlations always exist ex- 
cept for the cases (a) in which there is no dissipation (in 
agreement with the results of or (b) whe n the 

dissipation is maximal. For dissipative systems driven 
by thermal contact with a heat bath, it is not only the 
bath temperature which determines the steady-state of 
the system; the nature of the coupling to the bath is rel- 
evant as well (in contrast to non-dissipative systems, for 
which this last plays no role in the determination of the 
equilibrium state). With this in mind, we calculate the 
response of the system to a change in one of these defin- 
ing parameters, and compare the spatial dependence of 
the nonlocal response function to that of the correlation 
function. Finally, motivated by the question of energy 
flow in vibrated granular materials we study how 

the internal dissipation affects the nature of heat trans- 
fer, and calculate the energy along a chain whose ends 
are in contact with a heat bath. We note that the model 
studied in reference flSj in this context is the same as our 
model in the absence of dissipation. 

The system we study is a generalization of a model in- 
troduced in jl^ , which can be regarded as the mean- field 
version of the case of this Rapid Communication. We 
consider N particles localized on sites ti of a ID lattice, 
each characterized by its energy £"„. The entire chain is 
coupled to an external Boltzmann-distributed bath with 
temperature T. At each time step an interaction occurs 
in the system, either between system particles or with the 
external bath. Specifically, a particle is chosen randomly 
and its interaction follows the stochastic rule : 
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Here F is the overall rate of interaction of a particle; it 
sets the time scale and is irrelevant to the steady state. 
/, the strength of the coupling to the bath, is a constant 
which determines the probability of a particle interact- 
ing with the bath, a G [0, 1] is a parameter characterizing 
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the dissipation in an interaction (in analogy to a resti- 
tution coefficient), z is a stochastic variable uniformly 
distributed between and 1; and Eb is the energy of 
a particle chosen randomly from the bath. In what fol- 
lows, we shall only be interested in the steady states of 
the system. 

Our main results are: (i) correlations appear in this 
system for all < / < 1 provided a < 1; when f = 0,1 
or a = 1 there are no correlations; (ii) the spatial decay 
of the correlation function is calculated and found to be 
exponential in the limit of A'^ —f oo; (iii) the nonlocal 
response function (to a localized change in / or T) is 
proportional to the correlation function. This last point 
is reminiscent of the fluctuation-dissipation theorem of 
equilibrium statistical mechanics, but in our case there 
is no a priori reason to expect that the two functions will 
exhibit the same spatial dependence. 



CORRELATION AND RESPONSE 

Using the dynamics presented in |^ we can write equa- 
tions describing the time evolution of the moments of 
the energy distribution and the correlation functions in 
the system. We will be interested in the steady-state 
values of the moments, so we set the time derivatives 
to zero. It is simple to show that the average energy 
is given by En = E — 2-f-2a(i-f) ' '^ith T being the 



bath temperature in units where the Boltzmann constant 
is unity. In order to compute the correlation function 
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C{n,n + k) = EnEn+k — E , we need to calculate the 
second moments; this leads to a set of coupled equations: 
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1^ is a vector of constants. The N x N matrix A may 
be decomposed as A = T — B, where T is a tridiagonal 
matrix with constant diagonals (which can be inverted 
with the help of and S is an x matrix which is 
zero everywhere but the upper left 3x3 block (which is 
denoted by the 3x3 matrix Q). A~^ may be computed 
from the relation l2(l 



A-^ = T- 



(2) 



where Q is an x matrix of zeros except for the upper 
left 3x3 block which is the matrix {T — Q)^^ , X being 
the 3x3 unit matrix. In the limit A^ — > 00 this yields 
the result that the correlations decay as 



C{n,n + k) = D{aJ)e- 



(3) 



where D{a, /) is a continuous function, and A is the cor- 
relation length, given by = arccosh ^^^j^Z ■ 

We note that the correlation length diverges as a ^ 1 
and / — » 0, although we know that for a = \ there are no 
correlations; this is because D{a, /) ^ for these values. 
This means that around a = 1 and / = the correlations 
are the longest-ranged but the weakest. The behavior of 
the correlation length as a function of / and a is shown 
in Figure n in which results from simulating the model 
are presented for comparison. 
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In the last equation, k > 2. We note that if / = 0, 1 or 
a = I, then C (n, n + k) = 0: there are no correlations in 
these cases. This is consistent with the two-dimensional 
granular gas simulation of [l^ where velocity correlations 
disappear as the restitution coefficient goes to 1, and rem- 
iniscent of similar behavior of the density correlations [ll]. 
The case a = is unique in the sense that we obtain the 
mean field result for the distribution function and all cor- 
relations disappear except C(n, n + 1). 

For general / and a, the above coupled equations can 
be written as a matrix equation of the form AW = V , 
where is a vector whose fc*'' place is EnEn+k, and 



FIG. 1: The correlation length as a function of a for different 
values of /. 

Because of the correlations, the single-particle energy 
distribution function P{E) cannot be calculated ana- 
lytically for general a. For a = 1, detailed-balance 
holds, and an H-theorem may be proved, and the sys- 
tem comes to thermal equilibrium at the bath temper- 
ature T. When a = 0, the generating function g{uj) 
(defined by g{uj) = {e^'^^) = e-'^^P(E)dE) is the 
same as that of the mean field model 15]: 9{^) = 
{ujT + 1)2-Fi(l,2,2 - f,-ujT) where 2F1 is a hyperge- 
ometric function. 

For f — (and a ^ I) the system energy decays to zero 
and therefore P{E) = S{0) is the trivial steady-state. In 
Figure [3 we plot the energy distribution function P{E) 
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for different values of a. We note that the distributions in 
Figure pi are reminiscent of those found for the mean-field 
casehi 



pio 




FIG. 2: The distribution function on a logarithmic scale as 
found from Monte Carlo simulation with / = 0.5 and T = 1. 



is to a local change in temperature T at site n. As be- 
fore, we are able to invert the resulting matrix equation, 
which, in the limit N ^ oo |2Cl| gives: 
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In an equilibrium system, the response of an observable 
to a small change in its conjugate field is proportional to 
the correlations of this variable, with the system's tem- 
perature being the constant of proportionality. Although 
our response measurements do not take this form, it is 
interesting to note that in both cases, the response of the 
energy En+k to a local change in the bath interaction at 
site n is proportional to the correlation function Cn,n+k- 
Moreover, for given bath temperature, the ratios of the 
response to the correlation (for both types of response) 
decrease monotonically with a, and diverge as a — » 0. 



For a dissipative system coupled to a bath, both the 
coupling strength / and the bath temperature T deter- 
mine the behavior of the system. We will now consider 
the response of the system to a local change in the cou- 
pling to the bath. We stress that this is not a measure- 
ment which can be performed on a system in equilibrium, 
for which / plays no role in the steady state, and it is 
thus intrinsically a non-equilibrium measurement. Nev- 
ertheless, we shall see that, reminiscent of fluctuation- 
dissipation relations, the spatial response to such a per- 
turbation is proportional to the spatial correlation func- 
tion. 

We imagine that each site in the system is coupled to 
a Boltzmann-distributed bath at temperature T with a 
coupling strength /. We seek to calculate the change of 
energy at site n + k due to small change in / (say, from / 
to /o) at location n. Having changed the coupling at site 
n, the system no longer has spatial invariance, so instead 
of writing the time-evolution equations for a single par- 
ticle, we must consider the dynamics of the entire chain. 
The structure of the resulting eauation[2(ij| is again of the 
matrix form AE = V , where E is a vector of the ener- 
gies {Ej}, and V is a vector of constants related to the 
boundary conditions. A is a matrix which depends on 
a, f and fo which can be inverted to yield En+k, from 
which we obtain that, in the iV — > oo limit. 



(^)/o^/=i?(«,/)e-'=/^ 
«/o 



(4) 



where B{a, f) is a continuous function, and where A has 
the same value as in Equation This describes the non- 
local response of the system at a site a distance k away 
from the point of a local change in coupling to the bath. 
We note that this intrinsically non-equilibrium response 
has the same spatial decay as correlation function. 
We may similarly ask what the response at a site n + k 



HEAT TRANSFER 

We have thus far considered a system all of whose sites 
are in contact with the same external bath. As an appli- 
cation of the ideas of the previous section, we now turn 
our attention to the case where the ends of the lattice 
are coupled to heat baths at temperatures r_ and T+, re- 
spectively. For the case of a = 1, we expect that the tem- 
perature profile (i.e. the average energy at a site) will be 
a linear interpolation between the two wall temperatures: 
this is a manifestation of the Fourier law |l7l Il8l l2ll l22l| . 

The case where a 1 may be thought of as an ideal- 
ization of energy transport through an inhomogeneous 
driven granular gas, where we might imagine causing 
walls to vibrate and ask how this is manifested through- 
out the system. For example, Herbst, et al 5] and Barrat, 
et al 0] treat a collection of inelastic disks in 2D held be- 
tween vibrating walls, and find a characteristic decay of 
the granular temperature with distance from the wall. 
We modify our system so that of the particles on the lat- 
tice, only those at the end (fc = 0, -|- 1) are coupled to 
baths, at temperatures T_ and T_|_, respectively. For all 
other particles, the dynamics are as in Equation ^ with 
/ = 0. In similar fashion to the previous calculation, the 
relation for the average energy at site k may be expressed 
as a matrix equation, which may be inverted to give 
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T+sinh[{k -N - l)r?] + T_sinh[k7j] 
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(6) 



with T] = arccosh^^. This result is plotted in Figure 
13 We note that for a 1, the profile is linear, as 
expected^^. For all other a, sufficiently far from a wall 
the decay is exponential with decay length rj. It is inter- 
esting to note the similarity of these results to those of 
the numerical study of driven granular gasses in 
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This result may be understood in the context of the 
Fourier lawjl^l which relates the rate of heat flow, ^ 

to the local temperature gradient: ^ = —k^^- In our 
case, the energy Ei plays the role of both Q and T, and 
its time evolution is governed by: 
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The terms on the right hand side of this equation can 




FIG. 3: Average energy as a function of position for a 20 
particle lattice for different values of a. The symbols o are 
from a simulation of the model, the values from • from Equa- 
tionEl and the continuous curve is the solution of the damped 
diffusion equation obtained in the continuum limit. 

be interpreted as a sink term, and the heat flow out of, 
and into site i respectively. The net average energy flux 
into i is J- — —K{Ei — Ei-i) with heat conductivity 
coefficient k — 2{n+i) ■ note that in the continuum 
limit we obtain a damped diffusion equation for the 
steady-state: aV^E + 4(1 — a)E = with boundary 
condition E{x = 0) = aT+, E{x = N + 1)= aT-. 

In this Rapid Communication, we have described ex- 
act calculations for the spatial dependence of correlations 
and response for a model DDS. It is intriguing that, de- 
spite the fact that the fluctuation-dissipation theorem is 
not applicable to this system, the nonlocal response to a 
change in temperature or bath coupling has the same spa- 
tial dependence as the correlation function. Of course, it 
remains to be seen what the nature of the temporal be- 
havior of these functions is, and it would be surprising if 
they were to exhibit the same frequency dependence. 

We would like to thank Joseph Avron, Guy Bunin, 
Jean-Pierre Eckmann, Joshua Feinberg, Michael Fisher, 
Oded Kenneth, Yair Srebro, and Annette Zippelius for 
interesting and enlightening discussions. D. Levine grate- 
fully acknowledges funding from Israel Science Founda- 
tion grants 88/02 and 660/05, and the Technion Fund 



[9] 
[10] 

[11] 
[12] 

[13] 

[14] 

[15] 

[16] 
[17] 

[18] 

[19] 

[20] 
[21] 

[22] 



Granular Gases, edited by T. Poschel and S. Luding, 
(Springer, Berlin, 2001); Granular Gas Dynamics, edited 
by T. Poschel and N. Brilliantov, published as Lecture 
Notes m Physics 624 (2003). 

A. Kudrolh, M. Wolpert, and J.P. GoUub, Phys. Rev. 
Lett. 78, 1383 (1997). 

F. Rouyer and N. Menon, Phys. Rev. Lett. 85, 3676 
(2000); K. Feitosa and N. Menon, Phys. Rev. Lett. 88, 
198301 (2002). 

A. Barrat, E. Trizac, Phys. Rev. E 66, 051303 (2002). 
O. Herbst, P. Miiller, M. Otto and A. Zippelius, Phys. 
Rev. E 70, 051313 (2004). 

A. Puglisi, A. Baldassarri, and V. Loreto, Phys. Rev. E 
66, 061305 (2002). 

T.P.C van Noije and M.H Ernst Granular matter 1, 57- 
64 (1998); M.H. Ernst and R. Brito, Phys. Rev. E 65, 
040301(R); I. Goldhirsch, S.H. Noskowicz and O. Bar- 
Lev, Phys. Rev. Lett. 95, 068002 (2005); I. Goldhirsch 
and N. Sela, Phys. Rev. E 54, 4458 (1996). 

D. ben-Avraham, E. Ben-Naim, K.Lindenberg and 
A.Rosas, Phys. Rev. E 68, 050103(R) (2003); A. Santos 
and M.H. Ernst, Phys. Rev. E 68, 011305 (2003). 

E. Ben-Naim and P.L. Krapivsky, Phys. Rev. E 61, R5 
(2000); E. Ben-Naim and P.L. Krapivsky, Phys. Rev. E 
66, 011309 (2002). 

I Goldhirsch and G. Zanetti, Phys. Rev. Lett. 70, 
1619 (1993); M.R. Swift, M.Boamfa, S.J.Cornell and 
A.Maritan, Phys. Rev. Lett. 80 4410 (1998); S. Luding, 
R. Cafiero, H.J. Herrmann, Lecture Notes in Physics 624, 
293 (2003). 

D. R.M. Williams, E.G. MacKintosh, Phys. Rev. E 54, R9 
(1996). 

R. Soto, J. Piasecki and M. Mareschal, Phys. Rev. E 64 
031306 (2001); R.Soto and M.Mareschal, Phys. Rev. E 
63, 041303 (2001). 

A. Baldassarri, U. Marini Bettolo Marconi, and A. 
Puglisi, Europhys. Lett. 58 14 (2002). 

E. Ben-Naim and P.L. Krapivsky, Lecture Notes in 
Physics 624, 65 (2003). 

Y. Srebro and D. Levine, Phys. Rev. Lett. 93, 240601 
(2004). 

G. Peng and T. Ohta, Phys. Rev. E 58, 4737 (1998). 

R. Soto, M. Mareschal and D.Risso, Phys. Rev. Lett. 83 
5003 (1999). 

C. Kipnis , C. Marchioro and E. Presutti, J. Stat. Phys. 
27, 65 (1982). 

G.Y. Hu and R.F. O'GonneU, J. Phys. A 29 1511-1513 
(1996). 

D. Levanony and D. Levine, to be published. 

A.V. Savin, G.P. Tsironis and A.V. Zolotaryuk, Phys. 
Rev. Lett. 88, 154301 (2002); B. Hu, B. Li and H. Zhao, 
Phys. Rev. E 57, 2992 (1998). 

F. Bonetto, J.L. Lebowitz and L. Rey-Bellet, Mathemat- 
ical Physics 2000, Edited by A. Fokas, A. Grigoryan, T. 
Kibble and B. Zegarlinsky, Imprial College Press, 128- 
151 (2000). 



